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Abstract 

We study a general class of translation invariant quantum Markov evolutions for a particle on l"^ . We assume 
locality of the spatial jumps and exponentially fast relaxation in momentum space. It is shown that the particle 
position diffuses in the long time limit. We employ a fiber decomposition in momentum space of the evolution 
made possible by the translation invariance of the dynamics. A central limit theorem follows from perturbation 
. . . theory around a single fiber whose evolution in the momentum representation is described by a Markov jump 

process on the d-dimensional torus. 

o 

A classical problem in the study of dynamical systems is understanding diffusive behavior for 
^ some properly rescaled variable. When starting from a Hamiltonian dynamics, that often pro- 
ceeds in two steps. First there is a the identification of relevant space-time scales under which 
'"^ certain variables obey a reduced autonomous description. That specifies the limit starting from 
o I a microscopic dynamics and leading to a master or Boltzmarm-type equation, e.g. as the result 



1 Introduction 



of a weak coupling or a low density approximation, ||T[|2l,|3]]. Already there some irreversible 
^ behavior may be exhibited. Additionally, a second step can further specify the irreversible prop- 
g erties of a more restricted set of degrees of freedom. The present paper deals with the second 
i—i step, taking for granted a form of the master equation for the reduced description of a quantum 
^ particle hopping on the lattice. The specific derivation of that master equation, very much like 
>• a linear Boltzmann equation, starting from the unitary evolution of a particle in contact with a 
00 reservoir is not the subject of the paper, leaving background and further details to the literature, 
^ see e.g. [|4j|. 
(N ; 

Csj ■ We imagine a translation invariant law of motion wherein the free Hamiltonian evolution 
^ is interrupted by scattering events from interactions with the environment. The effective or 
resulting description is that of a Markovian open system. In quantum mechanics, Ltndblad 
equations take the place of Langevin or Fokker-Planck equations in classical probability theory 
describing a dynamical system under the influence of an idealized noisy environment, cf. |[T2|. 
This Lindblad equation is a aster equation for the evolution of the density matrix. The models 
' we study in the present paper are translation invariant Markovian evolutions for a quantum 
particle on . Its state is described via a density matrix pt for which, in position representation, 
the diagonal pt{x,x) gives the probability density for its spatial location x at time t. We show 
that under the right conditions it behaves diffusively, exactly the same result as for its classical 
counterpart: around the average position and as time t evolves, the spatial probability density 
Pt{x) becomes closer to a Gaussian, with parabolic space-time scaling a;^ ~ t after subtracting a 
possible systematic drift. We will also give some counterexamples to that result thereby making 
it less intuitive from a particular point of view. The questions and the applied techniques are 
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however quite similar to what has been studied starting from classical Boltzmann-type equa- 
tions, see e.g. (ZHSHH. A quantum example that is probably very related to ours is in IITOl . 

The next section gives a heuristic discussion of our model and of its origin, together with a 
first description of the main result. Section |3] repeats the model in a more general mathematical 
setting and the main theorem 13.31 is being stated. The following section discusses the essentials 
of the situation. The rest of the paper contains the proof of the main theorem and in a final 
Appendix some further explanation is given about the form and the proerties of the Lindblad 
generator. 



2 Description of the result 

We start directly from an open system dynamics, in terms of a quantum master equation. While 
we refer to ||4H51 for an introduction to the derivation of this quantum master equation, we give 
a heuristic description of the dynamics. We imagine the motion of our test particle to be steered 
by scattering events with a background fluid. The background fluid is homogeneous and the 
interaction with the particle is translation invariant. In a discrete time version we would say 
that the particle moves in the direction of its momentum where each time its energy and its 
momentum gets updated via the interaction with the background. In other words, first, if the 
momentum transfer is 9, we get the shift in momentum or boost 

p^e-^'^^pe^^^ (2.1) 

with X the position operator. In the rest of the paper, we set h = 1. Secondly, the state of 
the particle changes at fixed momentum. The result is given in momentum representation as a 
multiplication 

p{ki, k2) — ^ Me{ki, ^2) p{ki, ^2) 

with some positive kernel Mg{ki, k2), that we do not yet specify further. Our mathematical as- 
sumptions are basically formulated on the kernel Mg(ki, k2). These assumptions will express a 
certain locality in the jumps and they also allow a sufficiently smooth relaxation of the momen- 
tum. Specific expressions can be found in |[5l, and we come back to its spatial representation (for 
a particle hopping on the regular lattice Z*^) for the examples in Section |43l and also in Appendix 

m 

Finally we also add the free evolution with some Hamiltonian H that only depends on the par- 
ticle's momentum Vt. The continuous time version gives us now an evolution of the density 
matrix pt,t > 0, given by the master equation 

j^Pt = i[H, p,] + ^(p,) - i{vI/*(/), p,}, / = identity (2.2) 

where ^* : B{e{Z'^)) B{e{l/)) is the adjoint of a completely positive ma ^ : Bi{t^{Z'^)) 
Bi{e^{Z'^)) of the form 

^(p)= [ du{e) e-''^^Me{p)e'^^ , (2.3) 
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where u is a positive a-finite measure on the (i— dimensional torus T*^ and the Mg's are com- 
pletely positive maps that admit a Kraus decomposition of operators that are functions of see 
also Appendix lAl 

The position distribution is found on the diagonal, (x, x), which is itself not Markovian. The 
main and mathematical result of the paper is the identification of natural assumptions under 
which the measure fit, defined by 

X&y/tR—tv 

for a drift velocity v E and for an arbitrary region R C W^, converges in distribution to a 
Gaussian law. 

3 Model and results 

3.1 A translation-covariant semigroup 

We start with the kinematics. On the Hilbert space we define the translation operators 

Ty,y e U^, and the (vector- valued) position operator as 

(r,/)(x) = /(x + y), (X,/)(x) = x,f{x) for / G f (Z'^), xeZ" 

Here, and in what follows, the subscript j = 1, . . . ,d refers to the components in Z'^ or T*^. We 
will often consider the space H in its dual representation, i.e., as L^(T'^) where T'^ is identified 
with [— 7r,7r]'^. For g E L^(T'^), we define the vector of 'momentum' operators = (Vlj) as 
multiplication by A; G T'^, i.e., 

n,gik) = k,gik), ke[-7T,nf (3.1) 

Although the i7/s are well-defined as bounded operators, they do not satisfy [Xj, Clj] = — i, nor 
does ri generate the translations r^. 

The space of trace-class operators is denoted by A density matrix p G 

is positive and normalized, Tr[p] = 1. Through a Fourier transform, we can identify p as an 
operator on ^^(T'^) with a well defined integral kernel p{ki, k2) since p is trace-class. We refer 
to this representation as the momentum representation. It is useful to change variables and to 
write 

[p],(fc):=p(A;-|,A; + |), k,p (3.2) 

where we will think of [p]p as fibers of the density matrix p, indexed by p G T'^. The equation 
(|3.2|) defines a map 

[],:B,{e{Z'))^L\T''). 

We will study this fiber decomposition with more care in Lemma 14. 1[ In particular, our condi- 
tions will ensure that the function p ^ [p]p can be chosen in C^(T'^, L^(T'^)). The zero fiber refers 
to taking p = 0. 
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Now to the dynamics. The time-evolution is given by a norm-continuous completely posi- 
tive semigroup on Si(£^(Z'^)). As a result, by [20], $f satisfies 

l^.ip) = LiMp)) (3.3) 

where L can be written in the Lindblad form 

Lp = i[H,p] + ^{p) - l/2{^*{I),p}, p e i3i(£2(Z'^)) (3.4) 

with H a self-adjoint operator and ^ a completely positive map, acting on i3i(£^(Z'^)), and \E'* is 
its adjoint, acting on B{f{Z'^)). 

We want the semigroup to be translation covariant 

^tir*pT^) = T*^tip)r^, xeZ'^ 

This implies that one can choose H = H(il), a bounded function of the vector of momentum 
operators Q and, by [15], ^ is of the Holevo-form (|2.3|) further discussed in Appendix |Al A 
crucial consequence of this translation invariant form is the invariance of the fibers, i.e., the 
dynamics gives rise to an autonomous evolution on each fiber; see also below under (|4.5|) . 



3.2 Assumptions 

There are basically two sets of assumptions, one having to do with the spatial locality of the 
dynamics, and the other with the dissipativeness. They are closely related by Fourier-transform. 
We here formulate these conditions and we discuss them more closely in Section |4] after having 
stated the main result of the paper. 

The locality can first be formulated in terms of the completely positive maps Mg, the measure 
dz/(-) and the dispersion law H. 

Assumption 3.1. [Locality] We assume that the completely positive maps Mg* are defined by the kernels 

Me{ki, k2) as 

(Mep)(fci, k2) = Me{k^, k2) p{k,, k2) (3.5) 

where the function Mg{ki, ^2) is twice continuously differentiable in both ki, ^2 i^^d once continuously 

differentiable in 6. The function H is assumed to be twice continuously differentiable. 

Moreover we assume that dz/(^^) in (|2.3|) is a probability measure whose Fourier transform vanishes at 

infinity: 

I dz/(^)e^^^ — > 0, xeZ'^ (3.6) 

Jfd \x\/0O 

It is instructive to translate this Assumption 13.11 to real space such that its link to locality 
becomes obvious. In position representation 

*(p)(xi,X2) = ^ N{xi,yi,y2,X2)p{yi,y2): (3.7) 

yi,y2&''- 

for some function A^(-, - , - , ■). Translation covariance of \E' is equivalent to the property 

N{xi,yi,y2,X2) = N{xi + z,yi + z,y2 + z,X2 + z), for all2 G Z'^ 



4 



Our locality Assumption 13. II can be satisfiecflif 



Y] {{xi-yi)^ + {x2-y2f)\N{xi,yi,y2,X2)\ — > 
^ — ^ 0:1-2:2/^00 

This locality can be compared with asking a finite variance for the jumps in a random walk and 
it is supposed to exclude superdiffusive behavior. 

We now come to the dissipativeness. To formulate that conveniently, we introduce the 
Markov generator A on L^(T'^) 

(Af) (k) = [ du{9) r{k - 9) f{k -9)- [ du{9) r{k, k + 9)f{k) (3.8) 

with the transition rates r(k, k') defined by, 

r{k,k + 9) ■= Mg{k,k)>Q. (3.9) 

The measure diy{9)r{k^ k + 9) gives the probability per unit time of jumping to the state k + 9, 
conditioned on being in the state k. As in (|2.1|) the parameter 9 plays the role of momentum 
transfer. Remember that r(-, ■) is a C^-function by Assumption 13.11 One has < c||/||i for 

any / G C{T'^) and c := ||r(-, and hence A is the bounded generator of a contractive (and 
positive) semigroup on L^(T°'), see the connection with the Hille-Yosida theory in e.g. ||T9ll ). 
Dissipativeness here refers to ergodic properties. 

Assumption 3.2. [Dissipativeness] We assume 

1. A has a simple eigenvalue with eigenvector V G L^(T'^). 

2. The eigenvalue is separated from the rest of the spectrum by a gap Ba, 

bA ■■= - sup Re (spec(A) \ {0}) > 

The assumptions above guarantee that the semigroup generated by A, i.e., e*^, relaxes expo- 
nentially fast to the stationary distribution V. For future use, we let Yt stand for the Markov pro- 
cess on T*^ generated by A and started from V. Using standard techniques for Markov processes 
on compact spaces, constructive conditions are available for guaranteeing the Assumption 13.21 
in terms of r{k, k') and u. We say more on that in Appendix |Bl The basic idea is however that 
one verifies the second condition of Assumption |3.2| by requiring 

inf /" dv{9)r{k,k + 9)>{], (3.10) 

Using the L^(T°') space enables nice continuity properties of the fiber decomposition, to be dis- 
cussed under Lemma |4lTl 

While the above are natural ergodicity or gap-assumptions, in fact for our result we need less. 
In particular, the exponential relaxation is not strictly necessary. We will however not describe 
that in detail. 

*Note that there is a slight arbitrariness in choosing u and since only the measures du{9)M0 (ki , k2) are defined by specifying N. 
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The following standard construction will be useful in the statement of one of our results. 
Consider the scalar product 



(/, g)v := j dkV{k)f{k)g{k), /, g G C{T' 



and let H-p stand for the Hilbert space which is the completion of C(T'^) with respect to the scalar 
product (■, Define the quadratic form A-p by 

(/, Arg)v = j dkV{k){A*f ){k)g{k), /, g e CiJ') (3.11) 

where A* is the adjoint of A acting on L°°(T'^). 

The operator A appears naturally in a perturbation set-up around the zero fiber. In fact, the 
evolution on the zero fiber is the Markov process generated by A. 

We finally ask some properties that appear directly linked with the notion of diffusion. We 
have of course already that the particle must be sufficiently localized since it is described by a 
density matrix p E Bi{P{7/)), but we also ask that the first two moments are well-defined in 
the following way 

Xjp, XipXj, XiXjp arein ei(£2(Z'^)), fori, j = 1, . . . , (3.12) 

Products of bounded and unbounded operators such as in (|3.12[) are to be understood as kernels 
of sequilinear forms with densely defined domains. For example, Xjp E B{i'^{Z'^)) (which is 
implied since Bi C B)) means that the form b{f, g) := {Xjf, pg) restricted to (/, g) E Dom(Xj) x 
£'^{Z^) is bounded, and thus is extendable to I'^il/) x f' {!/■). Note that by the boundedness of 
the form g) and by the definition of the domain of the self-adjoint operator Xj, we have that 

pg E Dom(Xj), for any g E I'^ip), (3.13) 

so that the product Xjp makes sense. 

We write in general pt for the solution of the equation (|3.3[) with initial condition p^ E Bi. We 
will show that there is a f e M*^ such that 

V = lim - xpt{x,x) (3.14) 

That can be strengthened to a weak law of large numbers. One can obviously force ^' = by 
requiring some additional symmetries. Getting our results does however not depend on these 
extra requirements. In particular, equilibrium conditions such as detailed balance are mostly 
irrelevant for the diffusive behavior around the drift, except when one wants for example to 
relate the diffusion constant to the mobility in linear response relations. 

3.3 Result 

We define T*^^) and T^^^ as, respectively, a vector and d x d matrix of operators on L^{T'^), by 

{T^^^f){k) = i{VH){k)f{k)+ [ dv{e)m^P{k + e)f{k + e), (3.15) 
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where mf\k) = — iIm(ViMe(/c, k)) (and Vi and V2 are the gradients with respect to the first 
and second variables of Mg{ki, ^2)), and 

(T(2)/)(fc) = du{e) (^{mf^ - + e)f{k + 6)- mf\k)f{k)^ (3.16) 



where 



= I [(V? + Vl)Me{k, k))^^^^^ + (V? + Vl)Me{k, k)) . 

We let Pq be the projection corresponding to the eigenvalue of A, the Markov generator 
defined in (|3.8)) , and take S the reduced resolvent of A at the eigenvalue 0, i.e., the solution of 

^(0 - A) = {0 - A)S = 1 - Po, SPo = PqS = 0. 

Finally, recall that V is the eigenvector corresponding to the eigenvalue 0. Hence V{k)dk is the 
stationary probability measure on the torus. In general the projection Pq is non-orthogonal and 
has the form Pq = \V){ljd\, where 1^^ is the constant function on T'^ with value 1. We use the 
notation {g, f) := Jj^dk g{k)f{k) for the pairing between / G L^(T'^) and g G L°°(T'^). To the 
time-evolved density matrix pt, we associate the measures pt on MJ^, defined by 

Pt{R) ■■= Tr[l^^^,,(X) Pt], for a Borel set RcR'' (3.17) 

where l^/j+^j {X) is the spectral projection of the vector of position operators X on the set \/tR+ 
vt C M"^. 

Theorem 3.3. Assume Assumptions 13.11 and 13.21 and let the initial density matrix po ^ Bi{i'^{Z'^)) 
satisfy (|3.12[) . The limiting velocity exists and equals 

V := lim - xpt{x, x) = -1(1^^, T^^^P) (3.18) 

The measure ptr defined as in (|3.17|) with v as in (|3.18|) , converges, as t ^ 00, in distribution to a 
Gaussian with covariance matrix a = f3 + (3^, where (3^ is the transpose of the matrix j3, given by 

P ■= -{i^,^ (r(2) - T«5T(^))P). (3.19) 

where the operators and as defined above. 

The truncated second moments converge to the covariance matrix a, i.e. 

lim pt{x, x){xi - tv.j){xj - tvj) = cr{i,j) (3.20) 

Equation (|3.20|) tells us that the covariance matrix a can be interpreted as the diffusion tensor. 
In the case that T^^^ = iVH and T^'^'> = 0, it reduces to another matrix that below we call a. In the 
proposition below, the matrices a and a are "non-negative" in the sense of real valued vectors: 

yv G R'^{v,crv), {v,av) > 0. 
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Proposition 3.4. Consider the diffusion matrix a appearing in (I3.19P . 

1. a is non-negative. 

2. Define the (vector) function 

C ■.= VH - {VH,V) (3.21) 
on T^. The real-valued d x d— matrix a with entries 

= i^"c?tEp[(C(F*)).(C(V'o)), + (Cm)).(C(>^o)).], (3.22) 

is non-negative. (As in Section \3.2\ Yt is the stationary Markov process on T'^ with stationary 
density V and evolving according to the Markov evolution generated by A, cf. (|3.gP ). 

3. Assume that for all w G M.^, the function k t-^ {wX{k)) does not vanish identically (we have 
written (-, ■)for the scalar product on M.'^). That is, all components of the velocity VH can fluctuate. 
Assume in addition that A-p, defined as a quadratic form in (|3.11|) , extends to a bounded and sectorial 
operator on H-p. Sectoriality means that there zs < 7 < oo such that 

\{f,lm{Ap)fU < -7(/,Re(Ap)/)p, for all f e Hp. (3.23) 

Then the matrix a is strictly positive (i.e., it has strictly positive eigenvalues). 

4 Discussion 
4.1 Symmetries 

We can erase drift terms and simplify the diffusion matrix a by imposing further symmetries. 
Define the linear and antilinear maps U and V acting on / G £^(Z"') as 



(f//)(x) = /(-x) and {Vf){x) = f{x) 

By assuming the symmetry 

U<^t{p)U-' = MUpU-'), p G i3i(£2(Z^)) (4.1) 
(or on the generator L of the semigroup $t), we have no drift i.e. 

= V = lim - xptix, x) (4.2) 
tyoo t ^-^ 

X 

If rather than the symmetry U , we assume that the noise term ^ satisfies the symmetry 

V^(p)V'^ = ^(VpV-^), (4.3) 
then the drift depends on the stationary state V and the kinetic term H as 

V = lim - xpt (x, x) = {VH, V). 

X 



8 



When we have the symmetry (|4.3|) , then = and the diffusion matrix a simplifies to a 
matrix that can be written in the form 



a{t,j)= [ dk [ dz/(0)P(A;)([ViV2Me](„)(A;,A;) + [ViV2M,](,vo(A:,A;)) 

POO 

+ dtErmYtMciYo)), + (ciYtMayom, m 







where the rightmost expression is defined as in (|3.22[) , and the first expression on the right is 



a non-negative matrix. We sketch the reasoning for this as follows. Since Mg is a completely 
positive map on i3i(£^(Z'^)), for each 9, it follows that Me(A;i,A;2) is the integral kernel for a 
positive operator on ^^(T'^). This moreover implies that for any v E W^, {v, [ViV2M6)](A;i, k2)v) 
defines an integral kernel for a positive operator. The expression 

dk dv{e) V{k) [Vi V2M0] (fc, k)v) 

jd 

is thus non-negative, since the diagonal entries in the integral kernel for a positive operator are 
non-negative. Symmetrizing [ViV2Mg](^ij){k, k) in the above expression makes no difference 
in the evaluation of expressions {v, av) for t> G M'^ but ensures that we have a symmetric and 
real-valued diffusion matrix. 



4.2 Idea of proof: perturbation theory 

The main feature of our translation invariant models is that the evolution generated by (|2.2[) can 
be decomposed along the fibers (|3.2[) , i.e., one can write [Lp\p = Lp[p\p for some operators Lp 
and the fibers of the density matrix, [pt\p, obey the differential equation 

j^[pt]p = Lp[p]p, (4.5) 

The expression for Lp can be determined as a quadratic form through a trace formula: for F G 

L°°(T'^), 

'T,[L*[^vn^F{nyp'^'')p] = {F,Lp[p]p)= [ dkF{k)Lp[p]p{k). (4.6) 
A simple computation shows that Lp = ihp + \E'p with 

{hpf){k) := {H{k-^-)-H{k + ^-))f{k) (4.7) 

{^pfm := / dv{e){rp{k-e,k)f{k-e)-fp{k,k + e)f{k)), (4.8) 

where rp{k, k') = Mki_k{k — |, + |). Note that rp=o{k, k') was simply called r{k, k') in (|3.9|) . 
The notation ^p is consistent with the translation covariance of ^, which implies that [^p]p = 

Since only the fibers around p = determine the diffusive behavior, this suggests using a 
perturbation argument to capture the essential dynamical properties of those fibers. Under our 
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assumptions on A = Lq, it has a gap between the zero eigenvalue corresponding to the sta- 
tionary distribution V and the rest of the spectrum which has strictly negative real part. Hence 
small and sufficiently bounded perturbations of Lq also have a gap. 

The following lemma gives a condition on p such that the function p i-^ [p]p is twice differen- 
tiable. This assures the existence of the first two moments. 

Proposition 4.1. Assume that a density matrix p G Bi{l'^{l/)) satisfies 

Xjp, XipXj , X^Xjp are in Bi {f (Z"^) ) for i,j = 1, . . . ,d. (4.9) 

(as in (|3.12|) ), then the function T'^ — > L^{T'^) : p t-^ [p]p can be chosen to be twice continuously 
differ entiable. 

4.3 Examples 

We mention a certain subclass of models satisfying the conditions for (|3.3|) and another class of 
models that do not satisfy our conditions and that possibly show a non-diffusive (e.g. super- 
diffusive) behavior. 

We consider first the case when the completely positive maps in (|2.3|) reduce to M^p = p, 
such that ^ becomes 

^(p) = [ du{e) e''"" pe-''"" 
The map \l/ then acts as multiplication in the position representation 

^(p)(a;, y) = ^{x - y)p{x, y) 

where (p is the Fourier transform of the measure v. This is similar to the noise term for Gallis- 
Flemming dynamics | [13| for a quantum particle in M that has been frequently discussed in the 
decoherence literature | [Tni24| . The simplest case is when z/ is a multiple of Lebesgue measure, 
for which (p{x, y) = c5x,y and the particle changes position only by being carried through the 
kinetic term. 

Like in [13J, the particle can be thought of as in contact with a reservoir at infinite temperature 
and the stationary density V of the zero-fiber Markov process is the uniform distribution. Since 
Vi'V2Mg(ki, = 0, the diffusion constant a takes the familiar form of the second term in (|4.4|) . 
That is, the diffusion matrix reduces to a, as defined in Proposition 13.41 

Thinking about adding spatial jumps we arrive at models where the noise resembles a 
simple symmetric random walk. Yet, that easily breaks Assumption 13.21 Statement 2), see also 
(|3.10|) , and the model can become super-diffusive when the kinetic term is included. As an 
illustration, we consider a one-dimensional model. Let \E' have the form 

^'(p)(xi,X2) = ^i^i^yi^y2,X2) p{yi,y2), 

where N{xi, yi, y2, X2) is determined by a function r(x), x G Z'^, with a positive Fourier trans- 
form through the equation 

Nixi,yi,y2,X2) = r{xi - X2) x[\xi - yi\ = l]x[\x2 -Z/2I = 1] 
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where x\'] is the indicator (1 or 0). In that case, by the formal derivation in Appendix |Al the 
measure u{0) can be taken to be Lebesgue measure and the functions Mg{ki, ^2) are determined 
by 

Me{ki,k2) = 4 f(6') cos(A;i) cos(fc2) 

^(p)(A;i,A;2) = 4 [ d6 f (9) cos{ki + 6) cos{k2 + 9) p{ki + 9, k2 + 9) (4.10) 
Jt 

The smoothness of f{9) = Zln e'"^r(n) depends on the decay of c. Notice that the Markov 
process of the zero fiber is generated by 

{Af ){k) = 4 /" d9r{9) (cos(fc + 9Yfik + 9)- cos{kffik)). 
Jt 

That is describing a random walk on the torus T where an escape from the point k occurs with 
rate 4 cos^(A;) {Jj d9 f{9)) and the jump size 9 is independent of k and has probability density 

f{9) 
Ird9m 

Assumption 13 .21 is now violated in that, for A; = ±|, 

[ d9Me{k,k) =4{ [ d9r{9))cos\k) =0 
Jt Jt 

This implies that there are degenerate stationary distributions of the form 

V,{k) = \6{^ - A:) + (1 - \)6{ - I - A:) (4.11) 

for < A < 1. The process is thus slow to leave the regions around k = We conjecture that 
for an arbitrary smooth probability distribution V on T, e*^V converges in distribution to 7^1/2/ 
i.e., to (l4lT1l with A = i. 

When the kinetic term H is zero, we observe that still the diffusion constant as such keeps 
making mathematical sense as 

a = 2 [ dkVi{k){did2M){k,k) = % [ d9f{9) 
Jt ^ Jt 

When the kinetic term H is non-zero and H'{±.^) 7^ 0, we conjecture that the model exhibits 
a behavior closer to being ballistic with a width in position that grows on the order of t rather 
than t2. The basic idea is that when the particle attains a momentum close to A; = ±|, then it 
tends to stay in or around that momentum for long intervals with infinite mean residence times. 
We are not sure about the physical relevance or importance of these models. 

4.4 A classical analogue 

There is a sense in which our present quantum problem differs little from an analogous classical 
problem that starts from a linear Boltzmann equation. We explain that analogy here. Among 
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other things it relates our results to the recent interest and work on diffusive behavior in sys- 
tems of coupled oscillators where energy transport can be understood as a wave scattered by 
anharmonicities, ISJ. 

Consider a stochastic dynamics with state space 5 = Z'^ x T'^ such that the probability density 

Pt{x, k) = Tt{P){x, k) evolves as 

^Pt(x, k) = h{k){Pt{x + l,k)- Ptix - 1, k)) 

+ [ {T{x',k';x,k)Pt{x',k') -T{x,k;x',k')Ptix,k)), (4.12) 

for an initial Po{x, k) = P{x, k), where T(x, k; x', k') > is a transition matrix describing the rates 
of Poisson timed jumps from (x, k) to (x', k') and the symbol j^^, stands for jj^ dfc X^j-g^d- 
When T satisfies T{x + z, k; x' + z, k') = T{x, k; x', k'), then the master equation (|4.12|) describes 
a formally translation-invariant evolution. From a classical physics perspective the evolution 
(|4.12|l is still somewhat strange as it involves both jumps in position (x) and in momentum (k) for 
the single particle's state evolution. Still one can ask the classical questions about its diffusive 
behavior. We must then show that the centered position (xt — vt) / y/i converges in distribution to 
a Gaussian with non-degenerate covariance matrix. We show here how this can proceed along 
the very same lines as for the (indeed, more general) quantum case. 

By translation-invariance the dynamics does not "feel" the location of the particle and it 
follows that the "momentum" dynamics is Markovian, or, the marginal probability, Pt{k) = 
YlxeZ'i Pt{x-, k), for the k variable obeys an autonomous first order equation: 



d ~ 



Pt{k) = I dk'{f{k',k)Pt{k')-f{k,k')Pt{k)), (4.13) 



where T{k', k) = J2z + ^'i ^' This is analogous to the zero fiber of our decomposition 
for the quantum dynamics. Defining [Pt]p{k) = J^xeZ'^ e'^^-Pt(x, k), the dynamics (|4.12[) operates 
as [Tt]p[P] = [rt{P)]p for some semigroup [F^ : ^^(T'^) ^ L\T''). 

The closest thing to a joint probability density for position and momentum in the quantum 
case is the Wigner function. Here we go in the opposite direction and we define a "classical den- 
sity matrix" from a joint distribution function P{x, k). That is, for a momentum representation 
kernel, we formally define: 

p{k„k,) = Y: e-^'^-'^^P{x,'^^). (4.14) 

In this case, the positivity of p is lost, since it is not in general the kernel of an operator with 
positive spectrum. The invariant fibers again correspond top = ki — k2 in the momentum 
representation. The dynamics can thus be written in its fibers 

^\PtUk)=ismip)hik)[PtUk)+ [ dk'{fp{k',k)[PtUk')-fo{k,k')[PtUk)) (4.15) 
at jfd 
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where Tp{k\ k) = Yliz&d e'^^T(x, k'] x + z, k). One notices similarities with the fiber decomposi- 
tion and with (|47|) . 

5 Proofs 

We need a technical lemma to deal conveniently with trace-class operators. We write Bi = 
Bi{H), for some Hilbert space H; in the application to our result H = 

Lemma 5.1. Let C e Bi and let X be a self-adjoint operator. If XC G Bi , then 

\im-{e'^^ - I)C = XC (5.1) 

where the convergence is meant in the sense ofBi. 
Proof. Let 

C = ^Xn\fn){9n\, /„, fi-n G 7^, J^|A„|<00 (5.2) 

be the singular value decomposition of C G Si. Both families /„ and (?„ are an orthonormal set. 
By the comment (|3.13|) and Cgn = Xnfn, it follows that when A„ ^ 0, then /„ G DomX. 

If there are only finitely many terms in the singular value decomposition (|5.2[) , then the con- 
vergence (|5.1|l is guaranteed by an application of Stone's theorem for each n. When there are an 
infinite number of terms in (|5.2[) , an extra estimate is required to bound the tail of the sum. 

Defining the projection P„ = J2n=i Idn) {9n\, then we have the following bound 

1 ^ 1 

+ ||XC(J - Pn)\\, + ||-(e'^^ - I)C{I - P^)||i. (5.3) 

17 

Our strategy will be to pick, for each e > a number such that the last two terms on the right 
are bounded by e, for N > and arbitrary 7. Since the first term, with N = vanishes as 
7 \ by the reasoning above, the claim will follow. 

Since XC G Bi, we have \\XC{I — Pn)\\i ^ 0, as X 00. To bound the third term on the 
righthand-side, we note that for A = ^{e'^^ - I) and B = 2X, we have < |Ap < and 
hence also D*\A\'^D < D*\B\'^D for D = C(l — Pn)- Since ■ t-^ is an operator monotone 
function, we have D*\A\'^D < y^D*\B\'^D. Hence 

\\AC\\i < Tr ^C*\A\^C < Tr ^C*\B\^C = \\BC\\i. 

□ 

We continue with the proof of Proposition 14. II 
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[Proof of Proposition I4.1[ ] 

Step 1 

Assume the singular value decomposition for C, as in (|5.2|) . Then 

k ^ T,{k) :=5^A.(e^i^/„)(A;)(e-i^^0W (5-4) 

nSN 

is a L^-function (since it is the product of two L^-functions), which depends continuously on p 
(since p t-^ is a strongly continuous group on L^(T'^)). It is straightforward to verify that 
[C]p := Tp satisfies our definition of the fiber decomposition (|3.2|) . In other words, 

C{k,,k2):=Tk,-k,{'^^^) (5.5) 

is a kernel for the operator C. 
Step 2 

We show that, under assumption (|4.9|), [C]p is actually in Ci. If both C and XC, CX are in 
then 

^[{^.,C}]p=^[a- (5.6) 

because of (|5.1|) . The continuity of ^ [C]p in p, follows from (|5.6|) and from the general argument 
in Step 1, with C replaced by |{C, X}. 

Step 3 

The existence and continuity of the second derivative follows by repeating Step 2, since {C, X}X 
andX{C,X} arein-Bi. 

□ 

The following lemma lays out the standard perturbation theory [|T8| that we make use of. 

Lemma 5.2. Consider a family of bounded operators L.p on a Banach space for p G T'^. Assume that the 
g^g^ Lp are bounded and continuous as a functions ofp in some neighborhood ofO E T'^. Finally, assume 

that spec(Lo), the spectrum of Lq, contains as an isolated point, corresponding to a simple eigenvalue 
Then, for \p\ small enough, the operator Lp has a unique eigenvalue Dp such that 

Dp = {p, tr[T(^)Po]) + {p, tr[T(2)Po]p) - (p, tT[T^^^ ST^'^ Po]p) + o{p^), (5.7) 

where T^^^ and T^^^ are, respectively, a vector of operators and a self-transpose matrix of operators (i.e. 
T^-' = Tjf), defined by the expansion 

Lp-Lo = {p, T(^)) + {p, T^^^p) + o(p2), 

Po is the projection corresponding to the eigenvalue of Lq, and S is the reduced resolvent of Lq, at the 
eigenvalue 0, i.e. the solution of 

5(0 - Lo) = (0 - Lo)5 = 1 - Po, 5Po = Po5 = 0. 

Moreover, spec(Lp) \ {Dp} lies at a distance o{jp)from spec(Lo) \ {0}. 
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For our model the projection Pq has the form Pq = \V) {ljd\, where V is the stationary density 
for the Markov process of the zero fiber and lird is the indicator function over T'^. The first-order 
perturbation term T^^^ is given in (|3.15|) . The second-order perturbation T^'^'> is given by (|3.16|) 
and it only depends on the map ^, not on H. With the assumption of (|4.3|) , so that the second 
term on the righthand-side of (|3.15|) is zero, we have the following explicit expressions 

tr[r«Po] = i [ dkV{k){VH){k), (5.8) 
tr [T(2)Po] = - f dk [ du{9)V{k)mf\k), (5.9) 

J^d Jfd 

tr[r(^)ST(^)Po] = - / dt dkC{k){e'^XV){k), (5.10) 

with the function as defined in Proposition 13.41 We have used that on the range of S, the 
integral dte*^° is well defined and equal to 5". 

The expression (|5.9|) is non-positive since the matrices ViV2Me(A;, /c), are non-negative, as 
explained in Section 14.11 Seeing that the expression (|5.10|) is non-positive is a little more tricky, 
but it helps to rewrite the righthand-side as the following: 

dt dk{e''^Hmak)v{k)= dmvimxiXo)]- 

The evolution has been shifted to operate on the observables which can then be rewritten in 
terms of the expectation E^? of the Markov process Yt started from the stationary distribution V. 
We show the non-negativity of this expression in the proof of Proposition 13.41 

[Proof of Theorem 13.31 1 

To prove convergence in distribution for yUj, we show potntwise convergence of the charac- 
teristic functions. 

¥.,,(7):= / d/x,(x)e-^ = Tr[p,e^7f(^-^*)] = e-'v^'^^(lTr.,[e*^(p)]^) (5.11) 

where the third equality makes use of the fiber decomposition of our dynamics. In particular, 
we used the relation (cf. the proof of Lemma I^Tj) 

■^qx^^qx^ =[C] CeB,. (5.12) 

For a fixed 7 the limit involves only small fibers p (xt^K which suggests using a perturbation 
argument around the zero fiber. By Lemma (|5.2[) , for a small enough neighborhood of p = 0, say 
U c T'^, we have that 

[e*V]p = Ppe'""' + (1 - Ppy""' [pU := L, - D,Pp (5.13) 

and lle'^p || = 0(e^*^) as t /" 00, for some 6 > satisfying 6 — 6^ — > as p \ 0. The norm refers 
to the operator norm on B{L^{T'^)) and Ba is the gap of the operator A (Assumption |3.2[). 
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We now show that, in L^CT^), 



This follows by combining (|5.13|) , the relation 

1 



■5(7:0"7)p 



Dp = i{p, v) - ^{p, ap) + o(p^) 



(which follows from Lemma |5^ , and the fact that 

P. 



[p]p 



(5.14) 



(5.15) 



(5.16) 



where the first convergence is in B{L^{T''-)) and the second in V-(T'^). The first claim of (|5.16|) 
follows again from Lemma |5^ the second is a consequence of Proposition 14. 1[ We have shown 



(|5.14|) . Since norm convergence implies weak convergence, in particular [e*^(p)]^ integrated 



against the indicator l^d converges to the desired value. Hence converges in distribution. 

We now prove the convergence of the first and second moments (|3.19|) . By (|5.11|) , and the 
usual connection between moments and derivatives of the characteristic function, we have 
(whenever the righthand-side exists). 



7 ^ Pt{x, x){xi - tvi){xj - tvj) 



1, 

t dpidpj 



+ 7^(l,e*^1poU ^(l,e*^1poU ) (5.17) 
p=o opi p=oopj p=o 



Note that since the operator Lp has two continuous derivatives, the operators Pp,Vp (defined 

is bounded for p e U. One can 



oo. 



above) and the eigenvalue Dp do also. In particular 
see that, 

sup||(-^)e*^''|| =0(tV'*), t 
p(.u opidpj 

(and a similar bound for the first-order derivatives). By Lemma 14.11 the function p ^-^ [p]p is 
and we obtain 



dpidpj 



peu t/co 



0. 



Hence, in (|5.17|) we can replace e*^*' with Ppe^^^, for large t, and we see the convergence to a 
using the expansion (|5.15|) . 

□ 

[Proof of Proposition |3.4[ 1 



Proof of Statement 1) In the proof of Theorem 13.31 we showed the pointwise convergence of 
the characteristic function ip^^, i.e.. 



dfit{x)e 



1x7 



'|{7:0"7) 



tfoo 



(5.18) 
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(see e.g. (|5.11|) and (|5.14|) ). Suppose there were a 7 G M'^ such that (7, (17) < 0, then, for large 



enough t, \ip^^ (7) | > 1, which is impossible since ip^^ is the characteristic function of a probability 
measure. This proves the non-negativity of the diffusion matrix a. □ 

Proof of Statement 2) Now we consider the non-negativity of the matrix a. We show that a is a 
non-negative matrix by showing that the expression {w, aw) for w G M°' is always non-negative. 
Using an unsymmetrized form for a, we can rewrite our evaluation as 

POO 

{w,aw)= dtEp[/(Fi)/(>o)], (5.19) 
Jo 

where f{k) := {w, C,{k)) is real valued. Let the stationary Markov process Yt be extended to all 
negative values of t. Then G{t) = Ep[/(y't)/(Fo)] is an even function and so (|5.19|) is twice the 



value of the Fourier transform G{z) of G{t) at 2; = 0. We will show that G{z) is non-negative 
valued. Using Bochner's theorem we just need to check that G{t — s) defines a positive operator 
on L2(M). Let r] G ^^(M), then 



dtdsf]{t)G{t-s)r]{s) = ¥.r[\ / dt ii{t) f {Yt)Y] 

Jr 

where we have used the stationarity property E,p[f{Yt-s)f{Yo)] = E,p[f{Yt)f{Xs)]. □ 

Proof of Statement 3) The strict positivity of a is established as follows. Let w G M'^. By the 
assumption that the velocity fluctuates, g{k) := {w, C(^)) ^ T^v satisfies 

where Ijd G H-p, the identity function on T"^, is the 0-eigenvector of A-p. Note further that 

iw,aw) = {g,iAp)-^g)p, (5.20) 

since one can easily check that h := {Ap)^^g G Hp by using the fact that Ap is bounded and A 
has a gap. Assume that a is not strictly positive. Then there is a w G M'^ such that (recall that g 
and hence h depend onw) 

{g,iApr'g)p={Aph,h) = 0. 

In particular, this implies that 

{Re{Ap)h,h)p = 0. (5.21) 

Since — Re(Ap) is a positive operator (as follows from the fact that Ap is a Markov generator), 
we get Re{Ap)h = 0. By the sectoriality assumption, it follows that also lm{Ap)h = and hence 
Aph = 0. Indeed, for A > and v G Hp, 

\{{h + \v),lm{Ap){h + \v))p\ < -f{{h + \v),Re{Ap){h + \v))p = \'^{v,Re{Ap)v)p 

and hence the 0(A)-term has to vanish on the lefthand-side for all v. Since the zero-eigenvector 
of Ap is unique by assumption, we obtain h = cljd, c G M. This leads to a contradiction with the 
fact that h = {Ap)~^g and (1^^, "P) = 0. Hence a is strictly positive. □ 
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Note that (|5.19|) is related to the familiar central limit theorem for Markov processes 

lim ^ / dtf(Yt)^Af{0,a) 



where the convergence is in distribution, a = rft Ep[/(Ft)/(lo)], and the function / satisfies 

J^, dkf{k)V{k)dk = 0. 



APPENDIX 



A Holevo form 

The form of the Lindblad generator as a translation invariant particle dynamics has been ana- 
lyzed by Holevo, [16J. We give here more details concerning the relation of the Holevo form 
with the intuition of a spatially defined walk. We start with some formal calculations for defin- 
ing the maps Mg. 

We look back at p.7|) that we can also write as 

^(p)(si,X2) = Tti[N{I ^p^)]{xi,X2) 

for some positive operator G B{i'^{Z'^) £^(Z°')) which is partial trace class and satisfies 
A^(a;i + z,yi + z, 1/2 + z,X2 + z) = N{xi,yi, 1/2, X2); is the transpose of p in the position basis. 
Since is translation invariant, we can write N{xi,yi,y2, X2) = R{xi — X2, Xi — yi, X2 — 2/2) for 
some function R : Z^^ C. We can therefore define the map Fg as 

Fe{xi-yuX2-y2) = -7^ e-''^'^'~^''^N{x^,yi,y2,X2), 

where in the above we have held xi — yi and X2 — 2/2 fixed. Then FQ{ii,i2) defines a positive 
operator kernel on £^(Z'^), which we will assume can be written as Fo{ii, £2) = J2a fe,a{£i)fe,a{£2) 
for some functions fe,a on Z'^. Now looking at the action of the kernel Fe{xi — yi,X2 — 2/2) as 
determining a map on Bi{L?' {Z''-)) 

^ Fe{xi - yi, X2 - 1/2) p(z/i, I/2) = J]] ~ ^1) P^V^^ ^2) feA^-i - I/2) 

yi,y2 yi,y2 a 

which is a Kraus decomposition. The fe,a, fe,a act through convolution which means that they 
are multiplication operators in the momentum basis: 

X] feA^i) P(^i' ^2) feA^^) = Mg{ki, /cg) p{ki, k2) 
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and we get back 



^(p) = / dee''''Me{p)e-'''' (A.l) 

B Dissipativeness of the zero fiber classical Markov process 

We give some further discussion on the Assumption 13.21 While that brings us to the general 
theory of Markov jump processes on compact spaces which is not the subject of the paper, we 
do not give full proofs here but we sketch the central arguments. 
First of all, the conditions 

1. Continuity and strict positivity of r(fc, k + 6), 

2. Decay at infinity of the Fourier transform of the measure u, i.e. Assumption p.6|l 

appear sufficient for the irreducibility of the semigroup e*^, i.e. the property that for all g > 

0,9^0 

Snpp{e'^g} C Supp{^} => Snpp{g} = 

Since r{k, k + 9) is strictly positive, the support of the measure u determines the possible jump 
increments for the Markov process. Recall that the dynamical system of repeated shifts by a 
vector 9 = {9i, - ■ ■ ,9d) G T'^ is ergodic if and only if G X, the set of points satisfying 

ni9i + ■■■ + nd9d ^0 on T (B.l) 

for any combination of integers rii, ■ ■ ■ , G Z for which not all rij are zero, see [231. We now 
argue that the measure v has its support on J due to (|3.6|) , the decay of the Fourier transform 
ify. Consider the subset Ti^s,si,- ,sa) C for si, ■ ■ ■ , G Z, s G T consisting of all points in 6* G T"^ 
such that 

Si9i H h Sd9d = s on T 

Notice that (|3.6|) requires that, for m G Z tending to infinity, 

where formally dp{s) = L c/z/(6'i, ■ ■ ■ , 6^^). It follows that dp cannot have atoms and hence 

{lis, Si, -,3 a)) = 0- Thecomplementof the set J is a countable union of the subsets 7(o,ni,...,nj) over 
rii ■ ■ ■ Ud G Z. Hence, X has full measure and the Markov process must be irreducible. 

Next we exploit the inequality (|3.10|) . We write A = K + R where K and R correspond to 



the two terms on the right side of (|3.8|) . i? is a multiplication operator whose spectrum lies 
on the negative real line and by condition (|3.10|) we have that supspec(-R) < 0. Assumption 
(|3.6|) implies that the integral operator K is compact. Therefore, by Weyl's theorem, the essen- 
tial spectrum of A and R are equal and hence the eigenvalues above the continuous spectrum 
spec(i?) form a discrete set (see e.g. p. 101 of |[22]|). 

To end the argument, we must check that the only purely imaginary eigenvalue is at and 
that it is simple. This follows from a Perron-Frobenius argument, e.g. Thm. 13.6.12 in [14], which 
yields the desired result upon using that A is bounded and that it preserves C{T'^). 
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